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Invariance properties of physical systems govern their behavior: energy conservation in turbulence 
drives a wide distribution of energy among modes, observed in geophysical or astrophysical flows. 
In ideal hydrodynamics, the role of helicity conservation (correlation between velocity and its curl, 
measuring departures from mirror symmetry) remains unclear since it does not alter the energy 
spectrum. However, with solid body rotation, significant differences emerge between helical and non- 
helical flows. We first outline several results, like the energy and helicity spectral distribution and the 
breaking of strict universality for the individual spectra. Using massive numerical simulations, we 
'-^ , then show that small-scale structures and their intermittency properties differ according to whether 

helicity is present or not, in particular with respect to the emergence of Beltrami-core vortices (BCV) 
that are laminar helical vertical updrafts. These results point to the discovery of a small parameter 
t-|-H , besides the Rossby number; this could relate the problem of rotating helical turbulence to that of 

ryj ■ critical phenomena, through renormalization group and weak turbulence theory. This parameter 

O ' can be associated with the adimensionalized ratio of the energy to helicity flux to small scales, the 

" ^ , three-dimensional energy cascade being weak and self-similar. 

Ph. I. INTRODUCTION 

Turbulent flows, ubiquitous in nature, defy analysis due to their inherent complexity because of their nonlinearities 
leading to the strong coupling of a very large number of modes. Scaling laws beyond second-order, the order dealing 
with energy distribution in an incompressible fluid, are not amenable to simple dimensional analysis unless the flow is 
self-similar. Self-similar behavior in turbulent flows has been observed until now only at large scale, in the so-called 
I inverse cascades where excitations injected at an intermediate wavenumber reach lower modes as time elapses; the 
• ■ inverse cascade of energy in two-dimensional Navier-Stokes turbulence (see e.g., Kraichnan and Montgomery 1980; 
McWilliams 1984; Smith and Waleffe 1999; Boffetta et al. 2000; Tabeling 2002; Boffcta 2007), or the inverse cascade 
of magnetic potential in two dimensional and of magnetic helicity in three-dimensional magnctohydrodynamics are 
but three examples of this phenomenon. The link between inverse cascade, which emphasizes the emergence of order 
out of an otherwise chaotic flow as far as small-scale properties are concerned, and singularity of the underlying 
primitive equations is not clear. For two-dimensional Navier-Stokes turbulence (2D-NS hereafter) the regularity for 
all times can be attributed to the conservation of enstrophy, or squared vorticity (|^| 2 ) with u> = V x u the vorticity 
and u the velocity field; but regularity in the case of magnetohydrodynamics (MHD), in two or three dimensions, is 
an open problem. For 2D-NS as well as for surface quasi-gcostrophic turbulence for a rotating stably stratified layer 
(Picrrchumbert ct al. 1994), the self-similar behavior in the formation of eddies at scales larger than the energy input 
scale has been recently attributed, when examining the scaling properties of iso- vorticity lines, to conformal invariance, 
i.e. local scale invariance through transformations that preserve angles but not distances (Bernard et al. 2006, 2007). 
In that case, a link found to percolation theory allows for the analytical determination of scaling exponents such 
as the fractal dimension of vortex clusters, although the consequences for classical statistical measures of turbulent 
flows, e.g., through scaling laws for correlation and structure functions, has not been clarified yet. However, the 
three-dimensional case is known to be much more complex. In three dimensional homogeneous and isotropic fluid 
turbulence the flow is not scale invariant and the knowledge of one scaling exponent does not allow for the prediction 
of the exponents for all orders. The search for self-similar quantities in three-dimensional turbulence is a long-standing 
problem; it would relate its study with critical phenomena and the out-of-cquilibrium statistics of systems with a large 
number of modes, and it would allow the use of tools - such as the renormalization group (Ma and Mazenko 1975) - 
from quantum field theory, condensed matter and statistical mechanics to further our understanding of such flows. 

It is well-known that, under the influence of a strong external agent such as gravity, rotation or magnetic fields, 
several features emerge: on the one hand, waves are present in the flow (gravity, inertial, or Alfven waves respectively) 
and the interactions between nonlinear eddies and waves are not fully understood. On the other hand, the flow becomes 
anisotropic and in fact tends to (although never reaches) a near two-dimensional state, with thin layers in the case 
of stratification, or columnar (Taylor) vortices in the case of rotation. Since two-dimensional and three-dimensional 
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turbulent flows have vastly different dynamics, one may ask in what way are the three-dimensional dynamics of a flow 
under such circumstances altered, and in what way they may possibly be linked to the much simpler two-dimensional 
case, simpler be it only because regularity for all time is proven as mentioned before. 

In this context, we propose in this paper a brief review of rotating turbulent flows and specifically of the role that 
helicity (the correlation between velocity and vorticity) may play in such a dynamics. Rotating flows have been studied 
in the laboratory and observed in atmospheric flows as well as planets and stars (Greenspan 1968; Pedlosky 1986); a 
recent review of experiments performed on rotating turbulence, from 1975 to the present and both in the decay and 
in the statistically steady cases, can be found in van Bokhoven et al. (2009) (see Table 1 of their paper). Rotating 
flows have been investigated as well with theoretical tools and models of turbulence such as Large Eddy Simulations 
(Sagaut and Cambon 2008), and with direct numerical simulations (for recent studies of the specific effects of rotation 
on shear flows or on convective flows, see respectively Jacobitz et al. 2008 and Zhong et al. 2009). For example, in the 
atmosphere, close-range observations (2 km) with large resolution (50 m) during the VORTEX campaign (Verification 
of the Origin of Rotation in Tornadoes Experiments; Markovski et al. 1998) using Doppler On Wheels radars with 
deployed range of at most 120 s revealed many fine-scale structures, observed under the name of multi- vortex cores 
(Wurman et al. 1996): a quiet (laminar) eye, 200 m wide, surrounded by a ring of debris (implying turbulence) of 
similar width; moreover, at times, several concentric structures were observed as well, interpreted as different classes 
of debris with differing radar reflectivities. It has been known for a long time that multiple vortices are present in 
a tornado, for example from examination after the fact of the damages caused by the passage of the tornado; it is 
only recently that observational techniques have been able to quantify dynamical variables and localize such small- 
scale intense structures. In a recent study (Wurman 2002) it was shown for example that such multiple vortices are 
persistent and surround the central eye, with shear across them up to 100 ms _1 , and with vertical acceleration that 
can be five times that due to gravity, an observation in a turbulent flow that could be related to intermittency, as for 
example in the modeling of the formation of rain droplets (Shaw and Oncley 2001). 

We therefore proceed to the analysis of rotating flows with an emphasis on helicity and intermittency. Section [TT] 
explores the nonlinear dynamics of helical flows, Section [IIII discusses the structures that emerge in strongly rotating 
turbulent flows and their intermittency properties, and finally Section [IVI is the conclusion. 

II. NONLINEAR DYNAMICS OF HELICAL FLOWS 
A. The conservation of helicity 

One of the most important and useful principles of physics is that of conservation laws linked, through the the- 
orem of Emma Noether (1918), to invariance properties of the underlying equations; indeed, energy conservation, 
corresponding to invariance through translation in time, led Pauli (and later, Fermi) to hypothesize the existence of 
neutrinos. Similarly, the conservation of linear and of angular momentum is associated respectively with invariance of 
the dynamical equations under spatial translation and rotation. Angular momentum conservation is a key ingredient 
in the understanding of the relative motion of celestial objects, from accretion disks to planets, stars and galaxies. 

Conservation of energy is invoked when explaining the observation of the distribution of excitation among a wide 
range of scales in a turbulent flow: the nonlinear coupling due to advection, a convolution terms when transformed to 
Fourier space, leads to the feeding of modes at all the scales available to the system (Kolmogorov 1941; see also Frisch 
1995). In the two-dimensional case, the energy flows to large scales as hypothesized by Kraichnan (see Kraichnan and 
Montgomery 1980 for a review), where friction as is the case for the atmosphere of the Earth will stop the cascade 
from accumulating on the gravest mode. In three dimensions the cascade of energy is towards small scales and is 
arrested by dissipative processes. The fact that helicity {Hy = (u • u>), the correlation between the velocity and the 
vorticity) is also conserved was discovered much later than for the energy (Moreau 1961, Moffatt 1969). Note that 
in quantum mechanics, helicity is related to the relative direction of the particle's motion and its spin, and lack of 
mirror symmetry goes under the name of chirality. Helicity is not definite positive, unlike energy; it is a topological 
invariant, representing the degree of knottedness of vortex lines as well as the twisting of vortex lines, and it is a 
pseudo-scalar: its sign depends on the frame of reference, either right-handed or left-handed (the symmetry group 
related to its conservation is discussed in Yahalom 1995). Helicity in turbulent flows can lead to drag reduction, and 
to better mixing of chemical components in helical coherent structures (Duqucnne et al. 1993, Zimmerman 1996). 
Helicity has been measured in the atmosphere (Anthes 1982, Davies- Jones 1984, Rotunno 1984, Lilly 1986, Markovski 
et al. 1998, Lewellen and Lewellen 2007) and has been invoked to explain the long life-time of tornadoes and super-cell 
storms because of weakened non-linearities when it is strong. 

The lack of definite sign of hclicity has interesting consequences: one sign of helicity, say at small-scale, can be 
dissipated and yet represent a source of helicity at large-scale of the opposite sign. This renders the interpretation of 
helicity dynamics more complex. The invariance of Hy in the absence of kinematic viscosity v appears rather clearly 



The interplay between helicity and rotation in turbulent flows 



3 



when writing the Navier-Stokes equations in terms of the Lamb vector C = u x u> with V = p + |u| 2 /2 a modified 
pressure: 

-+uxu + 20xu=-VP + i/V 2 u + F; (1) 

at 

mass conservation reduces to V • u = assuming incomprcssibility for simplicity with a unit constant density and F 
is an external mechanical force that drives the turbulence, mimicking for example a convectivc input of energy. In the 
presence of an imposed solid body rotation with frequency f2, the rotation axis will be chosen to be in the z direction: 
f2 = fii; finally, note that V is modified by the centrifugal term in the presence of rotation. The global amount of 
helicity in the flow can be modulated through the forcing term (or in the absence of forcing, by the initial conditions); 
forcing (and/or initial conditions) can be taken with random phases with a prescribed Fourier spectrum, and with 
a prescribed amount of helicity through the coupling of two random vectors (see e.g., Pouquet and Patterson 1978). 
Alternatively, some large-scale order, emanating from an instability, can be prescribed around a given scale 27r/fco, e.g., 
through well-known flows like the Taylor-Green flow (TG hereafter) that mimics the laboratory experiments in fluids 
between two counter-rotating cylinders (see e.g. Monchaux et al. 2007) or Beltrami flows such as the ABC (Arnold 
1972). In the case of the TG flow, global helicity is zero because of symmetries although there are strong fluctuations 
of the local helicity density, whereas Beltrami flows have their velocity and vorticity parallel (or anti-parallel) globally. 
Moreover, Beltrami flows are known to be unstable (Arnold 1972; Kraichnan, 1973), and their dynamics is chaotic 
(Hcnon 1966), making them good candidates for the dynamo problem of generation of magnetic fields by turbulent 
flows (Arnold and Korkina 1983, Galloway and Frisch 1984; Gilbert 1991, Galloway and Proctor 1992) including at 
small magnetic Prandtl number (Mininni 2007). 

How helical a flow is can be evaluated with the relative helicity corresponding to the degree of alignment of the 
velocity and vorticity, i.e. i?y/(|u||a;|); it is defined in Fourier space as, with r(k) < 1 because of Schwarz inequality: 

r(k) = \H{k)\/kE{k) . (2) 

Turbulent fluids are known to develop helical structures which are persistent since their associated non-linear advection 
(the Lamb vector) is weak with u and us almost aligned (Pelz et al. 1985, Moffatt 1983, Moffatt and Tsinober 1992, 
Holm and Kerr 2002); thus, the temporal evolution of helical flow structures at the onset of the dissipative range where 
vorticity is strongest takes place on the slow dissipative time scale. Moreover, even when the total helicity is negligible, 
local helicity density is produced in the flow (Sanada 1993). This overall tendency toward alignment in the flow is 
fast (Matthaeus et al. 2008), occurring in a turn-over time tnl = Lq/Uq with Uq and Lq a characteristic velocity and 
length scale for the flow; it can be related to the alignment of pressure gradients and shear and has been observed both 
in DNS and in solar wind data. Even though small-scale structures are found to be strongly helical (of either sign), the 
helicity of a flow does not seem to alter its dynamics: indeed, evidence stemming from two-point closures of turbulence 
and from direct numerical simulations of isotropic, and homogeneous turbulence performed on incompressible fluids 
and with periodic boundary conditions, indicate (both for weak or strong global helicity in terms of r(k)) that the 
distribution of energy among scales follows a power-law (Kolmogorov 1941) which, expressed in terms of correlation 
functions or of structure functions of second order, reads ~ I 2 ' 3 , with S p (£) = ([ul(t + 1) — ul{v)] p ) the pth- 

order longitudinal structure function on a distance £, ul being the projection of the velocity field along the vector t. 
Similarly, it was shown using the rcnormalization group that the helical contribution to eddy viscosity is sub-dominant, 
with a r~ A dependence in the limit r — > oo, whereas the renormalized viscosity follows a classical r~ 2 diffusion law 
(Pouquet et al. 1978). 



B. The ideal case 



Kraichnan (1973) derived the statistical ensemble equilibria that emerge in three dimensions in the absence of 
viscosity for a system with a finite number of modes in terms of its global energy and helicity. Defining Fourier 
spectra such that Hy = (u • w) = J H(k)dk and Ey = (|u| 2 /2) = J E(k)dk, the equilibria are: 

47rfc 2 , , 87rfc 4 /? kd 

a(l — X) a A \\~x) a 

with a > and (3k max < a, k max being the maximum wavenumber of the truncated system; this condition is to ensure 
integrability, namely the positivity of energy, and \H(k)\ = 2\x\ = 2k\f3\/a < kE(k). The values of the Lagrange 
multipliers a and (3 are linked to the two invariants, Ey and Hy. In the non-helical case (/3 = 0), one recovers 
equipartition of energy among modes which, in 3D, leads to a k 2 energy spectrum. Equations ([3]) show that there is 
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TABLE I: Some of the possible types of cascades for fluid turbulence based on phenomenological arguments, taking into account 
the presence or absence of either helicity or rotation; Ttr is the transfer time and SI the imposed rotation. See text for details. 



J. Type / T tr — > 


TE, SI = 


TH, SI = 






K41 


tHcE 


D92 


Energy 


e=h=5/3 


2e+h=4 


e=h=2 


cascade 


a=2/3, b=0 


3a=7-3e, 3b=3e-5 


a=l/2=f=g, b=0 




c=-l/3, d=l 


3c=2(3e-4), d=-2b 


c=-l/2, d=l 


k D 


M/- 3 ] 1/4 




[tE/Sn] 1 ' 2 


Dual ? 


yes 


for e=h=4/3: K04 


yes 




B73 


K04 


M09 


Helicity 


e+2h=5 


e=h=4/3 


e+h=4, a+c=0, 


cascade 


3a=7-3e, 3b=3e-5 


a=l, b=-l/3 


3a+3b=2-f , b+d=l, 




6c=3e-7, 6d=ll-3e 


c=0, d=2/3 


3b=3e-5-2f, f+g=l 


k D 








Dual ? 


for e=h=5/3: K41 


yes 


for e=h=2: D92 



no accumulation of energy or helicity at small wavenumber, unlike the case of 2D-NS. The relative helicity grows with 
wavenumber although for k rnax — > oo and a > /3k max finite (corresponding to finite energy), one must have (3 — > 0. 

It has long been argued that inviscid dynamics is an indicator of what happens in the dissipative case. It is on the 
basis of ideal spectra peaking at large scale in 2D-NS that the inverse cascade of energy was postulated by Kraichnan. 
For three-dimensional Eulcr taking into account helicity, no such inverse cascade takes place in the non-rotating case. 
It was also shown recently that similarly to the non-helical case, inviscid dynamics leads to transient energy and 
helicity cascades that agree with the helical absolute equilibrium given above in Eq. ([3]): the excess of relative helicity 
found at small scales in the viscous run is interpreted as a thermalization (Krstulovic et al. 2009), and the large-scale 
spectra can be recovered with a model using a scale-dependent eddy viscosity. The ideal case in the presence of 
rotation is being studied presently. 

C. Dimensional scaling of Fourier spectra with or without rotation 

In the presence of viscosity, several spectral dynamics have been envisaged for helical flows. When generalizing 
them to include the effect of rotation, these possibilities are summarized in Table [H including some cases not derived 
previously in the literature. The various regimes are obtained assuming either (or both): 

• Constancy of the flux of energy £e = dEy/dt. 

• Constancy of the flux of helicity eh = dHy/dt. 

For simplicity, anisotropic effects are omitted (they can be added in a straightforward manner). The following 
timescales are also used, with te and th the turn-over times associated with energy and helicity (see below): 

• A characteristic time of the inertial waves in rotating cases as t w ~ I/O. 

• A cascade time of energy and helicity to small scales Tt r that can depend on a combination of te, th and t w . 

Note that the following type of arguments can be extended to other types of waves, such as for stratified flows, or 
in magnetohydrodynamics, for Alfven waves and whistler waves (see Galtier et al. 2005 for specific examples). 

Whereas the cascade time r tr is generally thought to be the eddy turn-over time te ~ [k 3 E(k)]^ 1 / 2 in the 
classical Kolmogorov cascade, one could envisage a characteristic time for transfer based on helical dynamics, 
th ~ [k 2 H(k)}^ 1 / 2 (Kurien et al. 2004), with te = th hi the case of maximal helicity, or alternatively a time 
that incorporates in some fashion the dynamical effects of waves. In that latter case, a phenomenological argument 
due to Iroshnikov (1963) and Kraichnan (1965) stipulated, in the context of Alfven waves in MHD, that the trans- 
fer time can be evaluated as r tr ~ te/c with e = t w /te the small parameter of the problem: in weak turbulence, 
the waves are assumed to be fast compared to the nonlinear coupling of eddies. This is in agreement with the fact 
that nonlinear dynamics is slowed down and weaker in the presence of waves; furthermore, the spectra evaluated in 
this manner coincide with those found in weak turbulence theory (Zakharov et al. 1992, Newell et al. 2001) when 
neglecting the effect of anisotropy (see Galtier et al. 2000 for MHD, and Galtier 2003 for rotation). 
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FIG. 1: Slices of vorticity in the xy plane in simulations of turbulence in periodic boxes. Top left: 512 3 simulation of non-helical 
non-rotating turbulence, Reynolds number Re = 1100. Top right: same with = 8 (Rossby number Ro = 0.07). Bottom left: 
1536 3 simulation of helical turbulence at early times with = 9 (Re = 5100 and Ro — 0.06). Bottom right: same at late times. 
Note the development of strong and smooth vorticity in the latter case, identified as columnar " Beltrami Core Vortices". In 
the top right plot, vortices associated to columns are barely observable, being less organized and with more small scale features. 



Specifically, what is assumed is that the energy and helicity spectra take the form: 

E{k) = C E e a E e b H nfk- e , H(k) = C H e c E e d H Wk- h , (4) 

with €e and 6h the energy and helicity flux rates respectively, Cg a generalized Kolmogorov constant and Ch the 
corresponding constant for the helicity spectrum, and with f2 the rotation frequency in the presence of an imposed 
solid body rotation as stated before. This choice covers the physics we examine in this paper and most cases studied 
in the literature, with energy, helicity and rotation all taken into account; the eight indices [a, 6, c, d, e, /, g 7 h] are the 
exponents to be determined through a combination of dimensional analysis and phenomenology, with e, h the spectral 
indices. 

The results of this analysis are summarized in Table U in its nomenclature, K41 stands for the classical Kolmogorov 
(1941) phenomenology, extended to a joint energy-helicity cascade; D92 stands for the energy cascade mediated by 
inertial waves (Dubrulle and Valdettaro 1992; Zhou 1995) and for its new extension to the helical case. B73 is for 
the dual cascade of energy and helicity as spelled out in Brissaud et al. (1973), including the case of zero energy 
flux. K04 stands for the case of a cascade fashioned by the helicity time-scale (Kurien et al. 2004). The case of a 
helicity cascade predominant over energy and mediated by rotation is M09 (Mininni and Pouquct 2009 ab). Finally, 
the generalization to the case where the time-scale based on helicity is the relevant feature of the cascade is dealt 
with in case tHcE without rotation (the similar case with rotation is omitted for simplicity but yields e + h — 3). 
In bold- face (K41, D92 and M09) are indicated cases that have been observed in the laboratory, in atmospheric 
flows or in direct numerical simulations; note that for a fc~ 4 / 3 helicity cascade (K04), the DNS observation concerns 
the so-called bottleneck effect between the classical Kolmogorov range and the dissipation range, an effect which may 
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FIG. 2: Left: Three dimensional rendering of vorticity using VAPOR (Clyne et al. 2007) showing the spatial juxtaposition 
of a Beltrami Core Vortex (BCV) and a vortex tangle, for the same run at late times as shown in Fig. 1 (bottom right). The 
width of the core is approximately 1/7 the size of the computational domain. Superimposed on the BCV are fluid particles 
trajectories which, in the vicinity of these structures, are laminar and helicoidal. Right: visualization at the same time of 
velocity- vorticity alignment. One can observe the core of the structure with positive (dark) relative helicity and concentric 
rings of alternate sign of helicity (dark and light) surrounding it, whereas, in the vortex tangle, no clear structure emerges. 



also be attributed to incomplete thermalization (Frisch et al. 2008) or prominence of non-local interactions near the 
cut-off wavenumber; fcjj is the dissipation wavenumbcr computed assuming equality of the viscous time l/[vk 2 ] and 
the relevant transfer time Tt r to small scales (either te, th or t^ l /t w = r^fl) at that wavenumber. 

The direction of cascades is not mentioned in Table [H it is well-known that the energy can be cascaded towards 
cither the small scales (standard three-dimensional case in the absence of rotation) or to the large scales (2D-NS case). 
It has also been observed to undergo both a direct and an inverse cascade in the presence of rotation (Smith et al. 
1996, Mininni et al. 2009, Mininni and Pouquet 2009a). Finally, by "dual" cascade is meant an energy cascade with 
a constant flux simultaneously with a helicity cascade with its own constant flux; this implies, using €e ~ kE(k)/rt r 
and en ~ kH(k)/T tr , the relationship 

H(k) = E(k) — 

£E 

which in turn implies e = h, f = g, c = a — 1 and d = b + 1 in Equ. ((4]) for this dual cascade. 

In more detail, the type of reasoning behind the estimations for the energy and helicity spectra listed above and 
detailed in Table [J is very simple: on the one hand, one argues that, by definition, an inertial range has constant 
flux, independent of wavenumbcr. A flux is the ratio of the total energy or helicity, divided by a characteristic time 
which we associate to the transfer time r ir . There are several candidates for r tr as discussed before. We eliminate 
the dissipation time which is supposed to be long compared to other relevant times in the problem, by definition of a 
high Reynolds number flow. 

So the transfer time can be the eddy turn-over time te ~ £/uf, this classical choice leads, in three dimensions, to 
the Kolmogorov energy spectrum and the helicity in this scenario follows the energy cascade with a constant flux as 
well. In the absence of rotation, that is actually the only solution that is observed: a dual Kolmogorov law (to within 
intcrmittency corrections), viz. 

E{k) = C E e 2 J 3 k~ 5/3 , H{k) = C H e E 1/3 e H k- 5 / 3 . 

This result of a dual Kolmogorov cascade was obtained using two point closures of turbulence (Andre and Lesieur 
1977), DNS (see e.g. Chen et al. 2003b and references therein) and the early phases of ideal dynamics (Krstulovic et al. 
2009). Note that, whereas the variation of Ce with Reynolds number has been documented with a slow convergence 
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FIG. 3: Zoom on vorticity intensity, when strong, at late times in the same simulation as in Fig. 2. The co-location of laminar 
structures, "Beltrami core vortex" (BCV) and of a tangle of vortex filaments with more complex paths and at smaller scale, to 
the left of the BCV, is clear. The BCV columns, different from Taylor columns (see text) are updrafts cyclonic vortices that 
are fully helical (see Fig. 2) and thus live for a long time, whereas the small-scale vortices are more evanescent, living for a 
couple of eddy turn-over times. Note two less-intense and inter- twined smaller multiple vertical vortex cores to the left. 



with Re (Ishihara et al. 2005), the corresponding variation of Ch with Reynolds number, or with different flows, 
is unknown at this stage and difficult to ascertain in the laboratory because of the difficulty to measure helicity. 
Furthermore, whereas for the ideal case, r(k) ~ k (see Eq. (|3]l). one finds numerically that in the dissipative case 
r(k) ~ 1/fc, i.e. indicating a recovery of mirror-symmetry in the small scales, albeit at a slow (1/fc) rate. Furthermore, 
although the maximal condition r(k) = 1 is not observed to be attained globally in turbulent flows, the alignment 
between velocity and vorticity is known to emerge in small-scale structures, namely in vortex filaments (Moffatt and 
Tsinober 1992, Matthaeus et al. 2008). 

Dropping the assumption of a constant helicity flux, one can find other solutions on dimensional grounds, namely 
c = 4/3 — h and d = h — 2/3. If instead of te a transfer time based on hclicity is used to regulate the cascade, 
other solutions can also be obtained. In the presence of helicity, a time that differs from te can be constructed using 
the sweeping time associated with the helicity spectrum, namely th ~ [^/(uetjJe)] 1 ^ 2 ', this leads to a fc -4 / 3 range, 
both for the helicity and for the energy (Kurien et al. 2004); in this latter case, one finds E(k) ~ e^e^^fe -4 / 3 and 

H(k) ~ e^fc" 4 / 3 (case K04 in Table BJ. 

As mentioned before, in the presence of waves a new characteristic time has to be taken into account as when 
t w < t e nonlinear interactions leading to a cascade of energy to small scales are damped. Weak turbulence theory 
(the small parameter of the problem that allows for closures of the equations being the ratio e = t w /te) leads to a set 
of integro-differential equations in terms of the various spectra of the problem. These complex equations can be shown 
to have both zero-flux (statistical equilibria) and constant flux (turbulent) solutions in terms of power laws of the 
wavenumber (for "warm' cascades that combine zero-flux and constant flux solutions for the energy, see Connaughton 
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FIG. 4: Left: perpendicular energy spectrum in a simulation of helical rotating turbulence with Re — 5100 and Ro = 0.06 
(solid), and in a simulation of non-helical rotating turbulence with Re = 1100 and Ro — 0.07 (dashed). The straight line 
indicates a power law ~ k~ 2 . Right: parallel energy spectrum for the same runs. In both figures the wavenumbers are 
normalized by the forcing wavenumber fcj?. 



and Nazarenko 2004). These weak turbulence solutions can in fact be recovered using a simple phenomenological 
argument which incorporates in a straightforward manner the small parameter e in the problem. Namely, one says 
that the transfer to small scales is slowed down as modeled by a longer transfer time Tt r ~ te/s- Taking t w ~ I/O 
for inertial waves, one obtains Tt r ~ t e /t w ~ Cl/[k 3 E(k)). Writing now that the transfer of energy to small scale is 
evaluated locally as €e ~ k 4 E 2 (k)/fl leads immediately to E(k) ~ [tE^Y^k" 2 . In order to obtain the dependency of 
the spectrum as written in Eq. ([][]), one uses dimensional analysis which yields a = f = 1/2 (Dubrulle and Valdettaro 
1992; Zhou 1995). The generalization to the case of a dual energy-helicity cascade is straightforward and gives for 
the helicity spectrum h = 2, g = 1/2, c = — 1/2, d=l as is given in Table [J (case D92); this results in: 

E(k) = C E [e E n] 1/2 k- 2 , H(k) = CH[n/e E } 1/2 e H k~ 2 ; 

the anisotropic version reads E(kj_,k\\) ~ [eE^k^^kj 5 ^ 2 and H(k±,k\\) ~ [eu / e]E(kj_, k\\), assuming that most of 

the transfer is in the perpendicular direction and that /cfj stems from initial conditions (one recovers isotropy with 

fcjj ~ fc|| ~ k±_). Note that this formulation of the spectra ensures that the dissipation wavenumber for the energy and 

helicity spectra is the same, namely [e^/i^fi] 1 / 2 (see Table IJ). 

What if one were to use a higher power of the small parameter e = t w /te in the evaluation of the transfer 
time? Let us take for example Tt r> 2 ~ te/c 2 . An analysis similar to what has been done above easily leads to 

E(k) ~ 4 /5 ft 4 / 5 fc- n / 5 , as proposed by Zcman (1994). The difference in terms of power law of the energy spectrum 
is small and given that intermittency may steepen the spectra, it will be hard to distinguish the Zeman spectrum 
from the D92 spectrum, viz. fc~ n / 5 versus k~ 2 . However, note that, by assumption, the dependency on the imposed 
rotation rate is significantly higher for the Zeman case, a point that may be amenable to examination using DNS and 
a parametric study; again, the generalization to a dual cascade gives h = —4/5 for the spectral index of helicity (note 
that this case is not incorporated in Table UJ). 

D. The helicity-dominated cascade to small scales 

One can make another type of hypothesis, supposing that, for some reason the energy transfer becomes negligible 
compared to the helicity transfer to small scales. This could be justified, e.g., on the basis that in the presence of 
strong rotation, energy flows to large scales and thus only a small amount of energy is available for a direct cascade 
to small scales. One is then faced solely with a helicity cascade, as done in Brissaud et al. 1973 (see case B73 in 
Tablc|T|. This leaves to an indctermination in the spectral indices of E(k) and H(k), with e+2h = 5, simply because 
the eddy turn-over time is expressed in terms of E{k) and the constant flux of helicity is in terms of H{k). This 
solution is compatible with the classical Kolmogorov case (e = h = 5/3) when both fluxes become comparable; it also 
admits e = 7/3, h = 4/3 which is the only solution with no dependence on ce but only on e#; it is also maximal 
(h = e — 1). A pure helical cascade ~ fc -7 / 3 is compatible with the exact law written in Gomez et al. (2000) based 
on the conservation of helicity when assuming maximal helicity leading to dimensional scaling (tu£ ~ ui/£ at scale £) 
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FIG. 5: Relative helicity in terms of the perpendicular (left) and parallel (right) wavemimbers in a 1536 3 simulation of helical 
rotating turbulence with Re = 5100 and Ro — 0.06. 



and assuming further no correlation between the velocity and vorticity fields. The difference of scaling in the presence 
of helicity may be related to the effect of large-scale flows (Olla 1998), as is observed for wall turbulence and the 
Bolgiano-Obukhov scaling for stratified flows; however, the lack of evidence for a —7/3 law in laboratory experiments 
and DNS is strengthened by the use of nearest-neighbor shell models of turbulence (Olla 1998). 

Although we are not aware of any observation relative to the e + 2h = 5 law put forward by Brissaud et al. (1973), 
the similar solution in the case of rotating flows has been observed in numerical simulations. The one important 
difference with the case of Brissaud et al. (1973) is that now the transfer of energy to small scales is mediated by 
waves and with a characteristic time taken to be r tr ~ t^/t w . This leads now to a new constraint on the energy and 
helicity spectral indices that read 

e + h = 4 , 

compatible with the dual cascade previously known (e = h = 2, sec D92 in Table [TJ but possibly leading to different 
spectra and thus to a loss of universality in helical rotating turbulence. The solution with e ^ h is clearly observed in 
several DNS (Mininni and Pouquct (2009 ab) and in Large-Eddy Simulations as well (Baerenzung et al. 2008, 2009). 
Why do we observe this interesting solution in the presence of rotation and not the corresponding e + 2h = 5 in its 
absence? The reason could be quite simple: at small Rossby number, the flow tends to a quasi two-dimensional state 
and thus supports an inverse energy cascade to large scales; the energy flux to small scale is thus negligible compared 
to the helicity flux to small scales, in accordance with the hypothesis of the models in this section. Note however that 
the energy that is transferred to small scales cannot be zero since E(k) > H(k)/k; evaluating this inequality at the 
smallest resolved scale of the flow, for example at the dissipation length kjj, indicates that the condition on the energy 
spectrum becomes smaller the higher the wavenumber available to the system, i.e. the higher the Reynolds number at 
a given rotation rate. Similarly, at a fixed Reynolds number, the higher rotation rate ensures more wave-dominated 
regime that leads to a more clear inverse energy cascade. This is consistent with results from a parametric study 
using a spectral model of small-scales that incorporates the effect of helicity on transport coefficients (Baerenzung et 
al. 2009). 

One can finally note that one may find a bit odd that, in a dual cascade, the sub-dominant field has in its expression 
a dependency on both 6e and e#. This is consistent, though, with the fact that finite dissipation must result for 
both spectra and that the dissipation is evaluated at ko whose expressions depend on the transfer time assumed to 
be relevant to the particular problem. 



E. A remark on polarized helical waves 

Several further comments concerning the analysis presented in Table U are in order. First, the constancy of the two 
(energy and helicity) cascades put sufficient constraints on the system for allowing a determination of most exponents, 
with some mild assumptions; otherwise, if only one cascade is hypothesized (say, the energy), a power law can occur 
in the spectrum of the other field (here, the helicity) that does not correspond to a constant flux solution for the 
helicity. Second, in some cases it may be difficult to define a power law solution for a non-definite positive field like 
the helicity, which may undergo changes of sign in its inertial range. A similar problem arises in MHD when looking 
at the cascade of the ideal invariant He = (u • b). Using second-order two-point closures of turbulence, it was shown 
that in fact the cross-helicity He, after a transient, establishes a two- lobe spectrum of one sign at large scale and 
of the opposite sign at scales smaller than the dissipation length (Grappin et al. 1982). A similar analysis has not 
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been performed for helicity in three-dimensional fluids and this remains an open question that can be tackled using 
both closure equations (Andre and Lesicur 1977) and DNS. However, when forcing the flow with Beltrami waves, 
the sign of helicity is predominantly that of the imposed forcing and thus one has no more difficulty evaluating the 
inertial index of the helicity spectrum in that case than for the energy. Furthermore, in MHD He can be seen as 
the difference of two (positive) energies, using the identity 4u • b = |u + b| 2 — |u — b| 2 , and its cascade can reflect 
subtle compensations between these two (pseudo-)energies |u ± b| 2 . Similarly, helicity can be expressed in terms of 
the difference of two energies. Indeed, the complexity of the behavior of helical flows can be better understood when 
considering simultaneously the [E,H] and the [i5 ± ,ff ± ] behavior of spectra where and H^ = kE^ refer to the 
cigenfunctions of the curl operator, corresponding to left-hand and right-hand circularly polarized maximally helical 
waves (Kraichnan 1973; Herring 1974; Waleffe 1993). Following Ditlevsen and Giuliani (2001) (see also Olla 1998), 
one can identify in the helical case a scale £r at which dissipation sets in by simply writing a balance between input 
and dissipation: 

ejj ~ v \ k 2 H(k)dk ~ vuui/£ H , 

where the upper bound in the integral is the dissipation wavenumber associated with the helicity spectrum, Kh = 
2tt/£h; using a dimensional estimate to ~ u/£ and the Kolmogorov scaling law u ~ e]^ 3 ^ 1 / 3 , one arrives at 

** - [- 3 4/4] 1/7 , 

a scale which is larger than the dissipative Kolmogorov length scale In ~ [e^/V 3 ] 1 / 4 ; this would imply that the 
helicity spectrum begins its dissipative range before the energy range in Fourier space. Note however that, if instead 
of supposing maximal helicity as done here, one simply expresses the dissipation in terms of the helicity spectrum 

H(k) ~ e E 1/3 e H k~ 5 / 3 

then one arrives at Kh = Kjj. The statement Kh/Kd < 1 is a point contradicted by numerical experiments 
(see e.g., Chen et al. 2003b, Alexakis et al. 2006). A remarkable explanation was given in Chen et al. (2003b) 
where one can find a detailed analysis of the interactions between all these modes based on the energy and helicity 
spectra as well as the spectra of the helical polarized waves: the scale £h determines the arresting of the spectra 
of the ± variables but the energy and helicity starts to dissipate at £o = 2n/k£>, the difference being attributed 
to cancellations between the fluxes of the ± waves. These considerations are related to the fact that, for maximal 
helicity, one would have H(k) ~ kE(k) ~ fc -2 / 3 , leading to divergences of its dissipation. It is a peculiar property 
of the helical wave decomposition that they are indeed fully helical (H^ = kE^) with H = H + — H~: helicity and 
its dissipation thus remain finite due to cancellations between the fluxes (Chen et al. 2003b). The fact that the 
dissipation scale for the helicity is equal to £jj is consistent with the idea of a balance between energy input and 
decay with H{k) ~ e^ 1 ^ 3 'enk~ 5 ^ 3 , as noted above. The scale £h enters in the helical wave decomposition to render 
the maximal helicity _ff ± constant. The dynamics of the [E, H] fields and the [£ ,± ,iJ ± ] fields thus differ, and the 
cancellations between the polarized waves occur in such a way that the helicity cascade is slaved to that of the energy. 
This remark also confirms the analysis in Kraichnan (1973) in terms of coupling of helical waves: a maximal helicity 
state is not consistent with the nonlinear dynamics even if in the initial state there is maximal helicity; in other words, 
a maximal helicity state assuming a Kolmogorov spectrum E(k) ~ k~ 5 / 3 leads to a non-physical helicity spectrum, 
which is only realized and observed in the ± variables. 

III. STRUCTURES AND INTERMITTENCY 
A. The emergence of Beltrami Core Vortices 

The preceding considerations rely on dimensional analysis on the classical basis of constant-flux solutions to the 
turbulence problem and using different time scales. Of all the solutions envisaged in Table [H only some have been 
observed, as stated before, and more work exploring fully parameter space remains to be done, with different Reynolds 
and Rossby numbers and different forcing functions. However, spectra are only a simple and constrained way to 
examine the data. What about the characteristic structures that emerge in rotating flows, and their statistical 
properties? 

We thus now report on some of the features of rotating turbulence obtained through a massive direct numerical 
simulation on a grid of 1536 3 w 3.6 x 10 9 points and forced with the ABC flow, with Lp = 2-K/kp the characteristic 
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scale of the forcing, of amplitude Fq: 



F 



Fq {[B cos(kpy) + C sm(kFz)] x+ 
+ [C cos(kpz) + Asm(kFx)] y + 
+ [Acos(kFx) + Bsin(kpy)] z} . 



(5) 



Results concerning the overall dynamics of the flow and its intermittency properties are reported in detail in Mininni 
and Pouquct (2009 ab). The Navier-Stokes equations in a rotating frame (see Eq. fT)) arc integrated with a pseudo- 
spectral code and periodic boundary conditions using a second-order Rungc-Kutta temporal scheme. The flow is first 
led to establish a statistically steady state with SI = 0.06, i.e. in the near absence of rotation, a phase taking roughly 
ten turn-over times; with Fq = 0.5 and kp = 7, the resulting rms velocity is U ~ 1 and the Reynolds number 



with the choice of v = 1.6 X 10 4 . Then, at a time labeled t = in the following, the rotation is set to f2 = 9, 
corresponding to a Rossby number 



The time step is At = 2.5 x 10 -4 . Note that the value of the forcing scale is chosen so as to let both a direct cascade 
to small scale and an inverse cascade to large scale develop simultaneously. In so doing, this run can be viewed as a 
combination of two computations performed earlier (Mininni and Pouquet 2009a) at lower resolution (with 512 3 grid 
points instead of 1536 3 here), in which either the direct cascade (with kp — 2) or the inverse cascade (with kp = 7) 
were studied separately. The computation is then performed for 30 turn-over times tjvl = Lp/U, corresponding to 
180 inertial wave periods. At early times, the dynamics is dominated by the Coriolis force and by resonant interactions 
between inertial waves. Because the dispersion relation is anisotropic and favors fcii = 0, the flow does not develop 
any substantial variation in the vertical and, in accordance with the dynamic Taylor-Proudman theorem, it displays 
columnar structures first clearly observed in the laboratory by Hopfinger et al. (1982) and numerically by several 
authors (see e.g. Bartello et al. 1994, Cambon and Scott 1999, Smith and Walcffc 1999). Note that the maximal 
vorticity in this flow is on average roughly 50 times its rms value and peaks are found up to 70 times that, whereas 
in a similar flow without rotation the maximal vorticity averages 20 times its rms value with peaks up to 30 times 
that value. The larger ratio of fluctuations to rms value in the vorticity of the run with rotation is associated to a 
substantial decrease in the enstrophy and the energy dissipation rates when rotation is strong. The maximum value 
of vorticity in the volume is however only weakly dependent on rotation, indicating that turbulence and small scale 
structures still develop in the flow even at small Rossby number provided the Reynolds number is large enough, and 
provided one waits long enough in terms of the rotation period. 

Figure 1 shows four horizontal slices of the vorticity intensity for a turbulent flow without helicity nor rotation (top 
left), a rotating flow without helicity at late times (t ~ 40, top right); a rotating flow with helicity at early times 
(t w 6.5, bottom left), and a rotating flow with helicity at late times (t s» 30, bottom right). The vorticity of a 
rotating flow without helicity at early times is not show as it looks similar to the non-rotating case. As a result, the 
vorticity in both helical and non-helical rotating flows at early times looks similar. As time evolves and anisotropy 
develops, the rotating flows tend to develop column-like structures in the velocity and vorticity (see e.g., Smith and 
Waleffe 1999). It is important to note, however, that such structures arc different in the helical and in the non-helical 
case. While the columns in the non-helical flow have small scale structure and display a myriad of vortex filaments, 
in the helical case a few columns develop a core with a smooth helical velocity field. One- of these structures can be 
observed on the left ( bottom- right slice) in Fig. 1. 

When visualizing the vorticity magnitude in three dimensions in a subset of the total volume for the helical rotating 
flow (see Fig. 2), using Perspective Volume Rendering as implemented with the VAPOR software (Clyne et al. 
2007), these laminar structures can be fully appreciated. Particle trajectories around the column are helicoidal, as 
materialized by the dark lines; the cores are fully helical and we thus name them Beltrami core vortices (BCV) to 
distinguish them from Taylor columns. These BCV are updrafts and cyclonic structures that live for long times 
(they were tracked for over ten turn-over times in the simulation) because their associated Lamb vector is weak; they 
support Kelvin waves along them, and move around because they arc embedded in the far-field velocity. Their number 
decreases with time through reconncctions because of the tendency of the flow to form larger-scale structures as time 
elapses in an inverse cascade. Overall, the vorticity is strongest in the tangle of vortex filaments, whereas the vertical 
velocity is coherent and strongest in the laminar columns (not shown). The vortex core of the column is surrounded 
by a calm region with weak vorticity (note the emptiness of vortical structures in the surroundings of the column) 
which acts as a transition region between the laminar and the turbulent flow. Far from these structures, column-like 
structures with a tangle of small scale vortex filaments also develop, as in the non-helical rotating case (see the dark 



R v = UL F /v « 5600 



Ro = U/{2VlLp) 0.06. 
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column on the left of the subvolume in Fig. 2). The BCV are Beltrami globally whereas beltramization in the tangle 
of vortices is local and random (see the relative helicity for the same region in Fig. 2). 

The origin of these stable structures can be identified by integrating backward in time the particle trajectories: they 
correspond to regions of large helicity where the columns form. Note that a rigorous scaling for the velocity profile 
in the coherent vortices that emerge at late times in 2D-NS is found to be u(r) ~ r~ 5 / 4 (Chertkov et al. 2009). How 
this scaling is altered by rotation is being presently investigated (Lebedev, private communication), and whether such 
scaling obtains in any of the coherent vortices that are found in DNS of rotating flows, either at early or at late times 
is an open topic for future research. The complex array of smaller-scales vortices next to the BCV is shown in more 
detail in Fig. 3, reminiscent of the vorticity field observed in many turbulent flows without rotation, except that one 
does perceive organization in the vertical direction in the form of a large-scale column, as is also found in the case 
of non-helical rotating fluids when the Reynolds number is large enough (Mininni et al. 2009). These columns have 
a more complex helicity structure and disappear on a time of the order of the eddy turn-over time. In the relative 
helicity of both structures shown in Fig. 2 (right) one can clearly distinguish the large-scale order associated with the 
BCV, with helicity of one sign, and a succession of rings of opposite signs (alternate light and dark regions), whereas 
as soon as one approaches the complex vortex tangle at smaller scale, no large-scale organization can be identified. 

The tangle of vortex filaments surrounding a laminar structure, together with more complex, larger and spiraling 
features, is reminiscent of observations of multiple core vortex tornadoes. Of course, only the very basic ingredients 
of helical convective storms and tornadoes are present in our computation: the input of energy that mimics the 
convective instability, the helicity often observed (Lilly 1986), and the strong rotation due to the local environment 
(Rotunno 1984). If obviously many other features of such an extreme event are absent in this bare-bone model (such 
as moisture, boundary effects, or micro-physical processes in general), there remains the possibility that the stability 
of these meteorological phenomena may be linked to the intrinsic dynamics of a rotating helical turbulent flow. The 
onset and ensuing acceleration of the tornadic motions can only be related to its environment and would require a 
more complete model; for example, it is already known that the growth of large-scale helical structures is observed in 
the presence of convection (Levina and Burylov 2006). 

The laminar organization of the velocity and vorticity fields is not observed in helical isotropic and homogeneous 
turbulence, nor is it observed in rotating flows without helicity as shown for example in Fig. 1. At early times, 
the bi-dimensionalization of the flow leads to columnar structures that are not surrounded by vortex tangles, the 
latter being due to the nonlinear terms coming into play. The interplay between rotation (which breaks the mirror- 
symmetry in the evolution equations) and helicity (which quantifies departures from mirror-symmetry of the flow) 
is the driving agent for the formation of strong localized and persistent columnar structures, even though helicity is 
mostly transferred towards smaller scales (Mininni and Pouquet 2009 ab) and is itself strongly intermittent. 

Finally, the next two figures show traditional statistics of turbulent flows, namely the energy spectrum (Fig. 4) 
reduced into its parallel and perpendicular components (i.e., in terms of variation of wavenumber parallel and per- 
pendicular to the rotation axis), and the relative helicity spectra (Fig. 5) also separated in parallel and perpendicular 
dependence. In the perpendicular energy spectrum the inverse cascade can be observed as energy piles up at wavenum- 
bers smaller than the forcing wavenumber. Noteworthy is the fact that a clear scaling at scales smaller than the forcing 
appears in the energy in terms of k±, with the non-helical case close to ~ kj 2 scaling (see e.g., Jacquin 1990, Zeman 
1994, Zhou 1995, and Cambon et al. 2004) and the helical case showing a slightly steeper spectrum (see Mininni and 
Pouquet 2009 ab). However, no clear scaling is observed in the parallel wavenumbers. In the helical run, an excess of 
relative helicity at small scales is observed in both the vertical and the horizontal directions (see Fig. 5). 

B. Intermittency and self-similarity in rotating flows 

Turbulent flows, in three dimensions and in the absence of rotation, display non-Gaussian statistics in their small- 
scales, associated with strong localized structures such as vortex filaments. How does the presence at small scales 
of laminar columns in rotating helical flows alter such properties? Specifically, intermittency is measured through 
scaling exponents of structure functions; for the velocity field, with ul the velocity component projected onto the 
distance r, one can examine: 

< 5u p L (r) >~ , 

assuming isotropy and homogeneity, and with the increment of a function / defined as Sf(r) = f(x + r) — f(x). In a 
rotating flow, one can also introduce increments on distances rj_ perpendicular to the axis of rotation and measure 



< 5u p L (r±) >~ r 
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FIG. 6: Scaling exponents f p in the direct cascade range of helical (stars) and non-helical (diamonds) rotating turbulence from 
direct numerical simulations. The dotted line corresponds to C, v = 0.71p and the dashed line to C, v = p/2. Note the deviation 
of the diamonds from the straight line for large values of p (see also Mininni and Pouquet 2009c). 



When ( p = a P, one talks of complete self-similarity of the flow with, in the case of the isotropic and homogeneous 
Kolmogorov (1941) energy spectrum, a = 1/3. Departures from such a linear scaling indicates that the flow is multi- 
fractal with a suite of exponents defining its small-scale properties and presumably characteristic of localized and 
strong intermittent structures; such departures can be computed exactly in the framework of a model of the passive 
scalar, but for three-dimensional Navier-Stokes turbulence, there is no such theory although some models arc quite 
close to the data, e.g. the She-Leveque model (1994). The observed curvature of the ( p = f(p) law is related to 
the fact that the probability distribution functions of velocity gradients have fat tails, with strong departure from 
Gaussianity at high values, corresponding to strong structures that are concentrated in space within the flow. No 
such departure from a linear law is observed for laminar flows, and similarly none obtains in inverse cascades toward 
large scales, as mentioned previously. 

In the presence of helicity, one may ask what are the intermittent properties of the flow in the absence of rotation. 
This was studied in Chen et al. (2003a) where they defined intermittency of helicity on its flux to small scales; they 
found that the intermittency of the velocity is comparable to that in the non-helical case, but that the intermittency 
of the helicity is larger than for the energy, insofar as the scaling exponents of high-order structure functions depart 
further from the linear scaling £ p = p/3, a result somewhat reminiscent of the case of the passive scalar. 

In the presence of rotation, several laboratory experiments and DNS have addressed this issue as well. New 
experiments using a non-intrusive (electromagnetic) forcing have allowed for detailed statistics of rotating flows using 
SPIV (Stereoscopic PIV) up to R\ ~ 240 and micro Rossby numbers of 0.15 (van Bokhoven et al. 2009). At low 
intensity of the forcing, self-similarity of the flow seems to obtain, with £ p ~ p/2, getting close to ( p ~ 3p/4 as 
the intensity of the forcing current is increased. Other experiments reported near £2 ~ 1 scaling (see e.g., Simand 
et al. 2000, and Baroud et al. 2002 for experiments with stirring), although small deviations from the self-similar 
scaling £ p ~ p/2 were also reported. As a result, whether rotation makes the flow self-similar or just decreases the 
intermittency is a matter of debate. In numerical simulations of forced rotating turbulence, £ p ~ p/2 with small 
intermittency departures (when compared with non-rotating turbulence) were found using non-helical forcing (Miillcr 
and Thiele 2007, Mininni et al. 2008). However, in the case of helical forcing at very large resolution, one obtains 
C p ~ 0.71p, and furthermore no deviations from the self-similar straight line was found within the error bars (Mininni 
and Pouquet 2009 ac). In the numerical simulations it was also found that the exponents get closer to the straight 
line as time evolves. Note that when self-similarity is observed, with £ p = ap, then a normalization of intermittency 
exponents by the data at some order n as in the ESS (extended self similarity hypothesis) methodology (Benzi et al. 
1993), gives trivially ( p /( n =p/n. 

For decaying flows, different intermittency exponents were observed as the flow decays (Morize et al. 2005, Sciwert 
et al. 2008), from classical non-rotating values at small times evolving toward £ p ~ p/2 as one approaches times of the 
order of l/fi, although strict self-similarity was not observed (the values found for the exponents in these experiments 
are consistent within error bars with the values reported in the numerical simulations of Muller and Thiele 2007, and 
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Mininni et al. 2008). For very long times in the experiments, £2 > 2, in contradiction with a E(k) ~ fc~ 2 law (D92 of 
Table HI, but in agreement with the law put forward in Mininni and Pouquct (2009a) that emphasizes the role helicity 
plays in the dynamics of a rotating flow; indeed, a strictly non-helical flow (implying orthogonality of the velocity and 
vorticity everywhere in space) is hard to obtain, even if, because of symmetries, the global helicity may be close to 
zero, as in the Taylor Green flow (Monchaux et al. 2007). 

Figure 6 illustrates these results showing the scaling exponents from two simulations of forced rotating turbulence 
with and without helicity. The exponents for the non-helical flow correspond to £ « 40 (fi = 8, and therefore 
tfl « 320). While £2 ~ 1, weak deviations from the self-similar scaling are observed for p < 5. For the helical flow the 
exponents were measured at t « 30 (CI — 8, and tCl w 270) and all exponents are consistent with ( p ~ p0.71 within 
error bars. 

Thus a possible explanation for energy spectra found to be steeper than k~ 2 in some of the experiments (e.g., 
Seiwert et al. 2008), and for the different behavior reported for the high order exponents, in some cases consistent 
with intcrmittency and in others with self-similarity, may be that some helicity is present in the flow, sufficiently so 
to affect the overall dynamics; this would be an interesting point to check both experimentally and numerically, in 
the latter case by varying the relative helicity of the forcing and observing the scaling that obtains. Another possible 
explanation may be related to some dependence of the exponents with the Rossby number, or with time. This point 
does raise the questions (see e.g., Seiwert et al. 2008, Baerenzung et al. 2009) of whether the transition from a 
Kolmogorov state to a steeper power law is smooth as the Rossby number is decreased, and whether for fixed Rossby 
number there is an asymptotic scaling the system reaches for very long times. 

In the case of helical rotating flows, the coexistence of laminar columnar Beltrami Core Vortices with a complex 
vorticity field at smaller scale must have an impact on the statistics of the velocity field; it is natural to try to associate 
the self-similar energy cascade illustrated in Fig. 6 to the dynamics of the BCVs whereas the helicity cascade (which 
is intermittent) may be associated with the small-scale vortex tangle; to that effect, a study using wavelets that allow 
to examine both the scale variation and the space variation of structures is in progress. 

The confirmation of self-similarity in turbulence in the combined presence of helicity and rotation (conditions that 
are relevant to many atmospheric flows) would relate the dynamics of such three-dimensional flows to the advances 
made in two-dimensional turbulence and critical phenomena in general. However, to use renormalization group 
techniques (Ma and Mazenko 1975, Forster et al. 1976) a small parameter needs to be identified, besides the Rossby 
number that governs the energetic exchanges between turbulent eddies and waves when dealing, e.g., with the weak 
turbulence regime (Zakharov et al. 1992, Connaughton et al. 2003, Nazarenko and Schekochihin 2009; see also Cambon 
et al. 2004, for a clarification of the link between two-point closures and weak turbulence). Indeed, the smallness of 
the ratio of the inertial wave period to the eddy turn-over time has already been used to derive integro-differcntial 
equations in terms of energy and helicity spectra (Galtier 2003) in the context of weak turbulence. However, the weak 
turbulence solutions are not observed in the numerical studies or in many atmospheric flows for at least three reasons: 
(?) the numerical resolution may be insufficient to see such laws; (ii) the theory is non- uniform in scale and the weak 
turbulence limit breaks down; (m) in the case of rotation the inverse cascade of energy is not present at lowest order 
in the theory and thus the solution selected by this approach is one of an energy cascade to small scales, whereas 
the numerical data of helical rotating turbulence indicates that this cascade is sub-dominant to the helicity cascade 
(Mininni and Pouquet 2009b). The candidate (backed by numerical data) for a small parameter in this problem is 
the (adimcnsionalizcd) ratio x = ^e/Lf^h, where 6e and e# are the direct energy and helicity fluxes, constant by 
definition in the inertial range. The energy flux to small scales is all the more negligible as more energy is transferred 
to large scales in an inverse cascade. Whether similar phenomena take place under the bi-dimensionalization of a flow 
due to other external constraints (such as stratification or an imposed magnetic field) is unknown at this point. 

IV. CONCLUSION 

The interplay between waves, nonlincaritics and energy /helicity transfer and intcrmittency is a topic of debate 
presently for a wide variety of waves (Connaughton et al. 2003) with, in some cases, the determination of power-law 
behavior in the wings of probability density functions (see, e.g. the discussion in the context of surface gravity waves 
in the ocean in Choi et al. 2005). However, the quasi-bi dimcnsionalization of the flow under an external agent such 
as rotation or a uniform magnetic field renders dimensional analysis more delicate since it introduces anisotropy in 
the scaling laws, and the presence of helicity has not been taken into account in general in these studies. 

The results presented here may indicate ways to study helical rotating turbulence from a theoretical point of view. 
In particular, we reviewed the many phcnomenological scaling laws derived in the literature and in the present paper, 
based on the different time scales present in the system, considering both rotation and helicity. Some of the scaling 
laws were reported in experiments or numerical simulations, while others where never observed. This may give us 
some information on what are the relevant time scales for the dynamics, as well as being a hint to what are the 
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dominant interactions between the waves and eddies. In particular, it is worth pointing out that in isotropic and 
homogeneous turbulence only the dual cascade of energy and helicity with Kolmogorov scaling has been observed so 
far (except for intermittency corrections and the bottleneck effect at the onset of the dissipation range), a solution 
which is obtained using the eddy turnover time and constant fluxes of both energy and helicity. In the rotating case 
more solutions have been reported, although all have the direct energy flux substantially decreased (and the time 
scale of the cascade increased) as a result of interactions with waves which select the resonances. 

Recent numerical simulations are consistent with different scaling laws in the direct cascade range of helical and 
non-helical rotating turbulence. We briefly compared these two cases and discussed the phenomenological arguments 
that are consistent with the observed spectra. Finally, we discussed recent experimental and numerical studies of 
intermittency in rotating flows. The conflicting results about intermittency (in some cases supporting intermittency, 
although decreased by rotation, and in others indicating self-similarity) may be associated to the effect of helicity in 
the flow, or to dependence of the intermittency exponents with the Rossby number or with time. However, when 
comparing numerical simulations, it was found that two runs at similar Rossby number and at similar times (albeit 
at different Reynolds number) display self-similar behavior or decreased intermittency depending on whether the flow 
had helicity or not. It is unclear for the moment whether the scaling exponents of the non-helical flow will behave as 
the helical ones for smaller values of the Rossby number or for later times. 

Self-similarity of these flows would open new possibilities of theoretical developments. It may lead the way to 
theoretical progress in unraveling the structure of turbulent flows, with possible extensions to the study of hairpin 
vortices in turbulent mixing layers and boundary layers (Rogers and Moin 1987) such as the planetary boundary layer, 
where stratification will play an important role as well. The finding of multi-scale structures that coexist spanning the 
range from the smallest dissipative scales in the flow to the largest energy containing scales, preserving scale invariance 
and with a small parameter associated with them, would relate the study of such complex flows to critical phenomena, 
where solvable models that preserve the complexity of the underlying processes exist. Finally, the self-similarity found 
at least for helical rotating flows can be exploited by subgrid models (as, e.g., in Baerenzung et al. 2009) which are 
often based only on second order statistics of the flow, to fruitfully study higher Reynolds numbers, lower Rossby 
number or larger separation of scales. Scale invariance provides the needed framework for the development of such 
models and should prove particularly fruitful in such cases. Several sub-grid models can be devised in this context. 
Lautenschlager et al. (1988) proposed, on the basis of similar analyses performed in the case of coupling to a magnetic 
field, to add, as parametrization of the small scales, an expression of the form au) +7AU;, with a and 7 depending on 
the amount of helicity in the flow, similar to the alpha effect in MHD whereby a large-scale magnetic field is unstable 
due to the small-scale helicity of the flow. It was shown in Pouquet et al. (1978) using the renormalization group 
that in fact the a term above is missing when performing a systematic expansion in terms of elimination of the small 
scales and that, indeed a term proportional to k 3 was present, although it was deemed negligible in the limit k — > 
when compared to the renormalization of the viscosity. In related studies, Frisch et al. (1988) showed that the AKA 
(anisotropic kinetic alpha) effect can be important for flows lacking parity invariance (invariance under simultaneous 
reversal of position and velocity), leading to destabilization of the large scales. This may be compared with the work 
of Yokoi and Yoshizawa (1993) for effects in inhomogeneous flows that can arise in the laboratory through imposing 
a swirling flow to a pipe flow, or in the atmosphere by combining swirling motions with updrafts linked to either 
convective motions or boundary layer effects. 

Finally, we believe several of the recent results suggest that a better understanding of the role of the inverse cascade 
of energy in rotating flows is required, and is a likely next step in the investigation. How does it coexist with the direct 
cascade of energy in the non-helical case, and the direct cascades of energy and helicity in the helical case? How does it 
compare with the two-dimensional case in the absence of rotation, for example in the statistical properties of vortices 
(see e.g., Macwilliams 1984)? The similarities and differences with the inverse cascade of energy in 2D-NS may help 
us understand the self-similar behavior, and relate the problem with the recent advances in conformal invariance and 
two dimensional turbulence (see Bernard et al. 2006 and Cardy et al. 2008). 

Computer time provided by NCAR which is sponsored by NSF. PDM is a member of the Carrera del Investigador 
Cientifico of CONICET. 
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